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ABSTRACT 


L characteristic length 


As part of its development of cryogenic 
fluid management techniques for spacecraft, 
the NASA Lewis Research Center Cryogenic Fluid 
Technology Office (CFTO) is planning to 
perform ground tests of nonvented fill tech- 
niques on a 4.96 cubic meter light weight 
liquid hydrogen tank. This tank is similar in 
size and shape to the tankage planned for 
CFTO 1 s COLD-SAT liquid hydrogen flight ex- 
periment. This paper presents the analyses 
used to select two injection systems for 
nonvented fills of this tank at design flow 
rates between 220 and 450 kg/hr. The first 
system uses multiple nozzles spraying from the 
top of the tank through the ullage space. 
This system should be capable of liquid fill 
levels in excess of 95 percent. The second 
system injects the liquid through a submerged 
nozzle and should produce fill levels on the 
order of 80 percent liquid. 


M mass 

M mass flow rate 

n number of drops in spray 

Nu Nusselt Number 

R radius 

t flight time 

Re Reynolds Number 

u internal energy 

v velocity 

V volume 

w rate of work 


NOMENCLATURE 

a tank semi-major axis 
A area 

b tank semi-minor axis 
C 0 discharge coefficient 
C p specific heat 
D diameter 
G volumetric flow rate 
h enthalpy 

h heat transfer coefficient 
k thermal conductivity 


x liquid fill height 
a thermal diffusivity 

X scale factor 
\i viscosity 
p density 
a surface tension 


Subscripts 


cond 


condensation 


fullscale prototype OTV value 


gas 

H2 


gaseous property 
hydrogen property 


1 


in 

inflow property 

inf 

interface 

ig 

between liquid and gas 

liq 

liquid property 

model 

value for tank being scaled to (RPM 
tank) 

sat 

saturation property 

sgas 

saturated gas property 

system 

system property 

wall 

tank wall property 

water 

water property 

INTRODUCTION 


On orbit transfer of cryogenic liquids is 
considered enabling to many future NASA mis- 
sions, from space based orbital transfer 
vehicles (OTV) to manned mars exploration. 
The techniques required to transfer cryogens 
in low gravity are quite different from those 
used terrestrially. During a normal gravity 
fill a top vent is kept open to vent the vapor 
generated during the fill process thereby 
maintaining a low tank pressure. If the 
normal gravity technique is used on orbit, the 
uncertainty of liquid and vapor distributions 
in low gravity may result in dumping of large 
amounts of liquid overboard. The No-Vent Fill 
process is a methodology used to reduce fluid 
loss by allowing the tank vents to be kept 
closed while the tank is filling (Chato, 
1988) . 

The procedure works as follows: 

The tank wall is cooled to a temperature 
sufficient to remove most of the thermal 
energy from it. The tank is filled with a 
subcooled liquid so that the end state condi- 
tion is also subcooled liquid. Sprays and 
Jets are used to mix the fluid and maintain 
conditions in the tank close to thermodynamic 
equilibrium. Providing the initial energy in 
the tanks is low enough, the equilibrium 
mixture keeps the tank pressure low so that 
the entire tank can be filled without opening 
the vent valve. 

The concept of No-Vent Fill has been 
extensively analyzed (Merino et al, 1978; 
Wiilen et al, 1981), but very little testing 
has been done (see Fester et al, 1970 for 
liquid fluorine testing) and, to the authors 
knowledge, no liquid hydrogen data has been 
published. To obtain a more empirical under- 
standing of the No-Vent Fill problem a series 
of liquid hydrogen experiments is to be con- 
ducted at the NASA Lewis, Plum Brook K-Site 
cryogenic vacuum chamber facility using an 
existing lightweight liquid hydrogen tank. 
This paper presents the rationale used to 
design the fill systems for the existing tank- 
age and their predicted performance during the 
No-Vent Fill process. 


SUMMARY OF THE EQUATIONS FOR THE NO-VENT FILL 
PROCESS 


In his previous work (Chato, 1988) the 
author developed a model of the No-Vent Fill 
process. This model separates the fill pro- 
cess into two stages. The first stage is the 
"Liquid Flashing Stage". This stage is chara- 
cterized by flashing and boiling and is model- 
ed as an equilibrium energy balance between 
the hot tank wall and the incoming liquid 
flow. The equations for this phase are: 


Wall energy balance (assuming the liquid 
inflow vaporizes on striking the wall but does 
not superheat) 


= M in (h sgas - h,„) 
dt 

Gas mass balance 

(1) 

dM = M 

^^gas in 

dt 

Gas energy balance 

(2) 

^gas^igas ~ ^in ^sgas - U gas) 

dt 

( 3 ) 

The second stage the "Vapor Condensation 
and Compression Stage" divides the tank into 
three nodes; gas, liquid, and interface. It 
then models the energy transport between them. 
The equations for this stage is as follows: 

Gas mass balance 


^^gas ” ”^cond 

dt 

(4) 

Gas energy balance (neglecting gas phase 
heat transfer) 

M gas dUg as + M C ond U gas ^cond^gas + W lg 

dt 

( 5 ) 

Liquid mass balance 


M cond 

dt 

(6) 

Liquid energy balance 


^liq— — liq + U liq^^liq + Qjrtf + ^in^in + ^cond^lg 

dt dt 

II 

K • 

IQ* 

Interface mass balance (assuming an 
thin interface) 

infinitely 

^Minf = 0 
dt 

Heat transport from bulk 
interface 

(8) 

liquid to 

Qinf = hA inf( T sat” T liq ) 

( 9 ) 

Interface energy balance 


^cond — — 1 iq-i- 

< h 3„s- h li q ) 

(10) 


2 



(13) 


Compression work 


For the RPM tank: 


w lg = _Pg as (M (n + M cond ) (11) 

P 1 iq 

These equations are solved with a com- 
puter algorithm called NVFILL which uses a 
finite difference approximation of these 
equations. Refinements have been made to the 
NVFILL algorithm to model the test tank which 
will be discussed below. 

EXPERIMENTAL HARDWARE 


X = 0.489 

Heat transfer correlations for the No- 
Vent Fill processes are not well developed, 
but a couple of scaling estimates can be made 
based on the assumed behavior of E, the inter- 
facial heat transfer coefficient. Defelice 
and Aydelott (1987) assuming a constant heat 
flux per unit area (and hence constant E with 
the same driving temperatures) give the scal- 
ing criteria for inflow as: 


Test Tank 


^model 


‘full scale 


(14) 


The test tank selected was one designed and 
built for use in the Research Propulsion 
Module (RPM) program conducted by Lewis Re- 
search Center in the early 70* s (DeWitt and 
Boyle, 1977). The RPM liquid hydrogen tank is 
ellipsoidal with a 222.5 cm major diameter and 
a 1.2 to 1 major to minor axis ratio. The two 
ends are joined by a short 3.81 cm cylinder 
section. The tank is made of 2219 aluminum 
chemical milled to a nominal thickness of 2.21 
mm. Thickened sections exist where required 
for manufacturability (mainly weld lands) . 
The tank has a 71 cm access flange on the top. 
Tank weight is 149.35 kg. Tank volume is 4.96 
m 3 . The tank is covered with a blanket of 34 
layers of multi-layer insulation (double 
aluminized mylar with silk net spacers) and is 
supported by 12 fiberglass epoxy struts. 
Thermal performance of the tank is documented 
in DeWitt and Boyle, 1977. Figure 1 shows the 
tank installed in its support structure. This 
tank has several features which make it desir- 
able as a test bed for spacecraft technology: 
It is of the same lightweight chemical milled 
construction used in space flight tanks. It 
has a multi-layer insulation blanket with 
performance nearly identical to current in- 
sulation designs for Orbital Transfer Vehicles 
(OTV) . In addition, the tank is similar in 
size and shape to the tankage planned for 
CFTO's COLD-SAT liquid hydrogen flight experi- 
ment. This will allow the CFTO to asses many 
of the problems expected to be encountered in 
the COLD-SAT spacecraft, even though many 
technologies will require flight test in low- 
gravity for their ultimate proof of concept. 


This scaling is probably most applicable to 
highly mixed cases where the heat transfer 
coefficient asymptotically approaches some 
maximum value. As an alternate bound for 
cases where a large amount of stratification 
is evident an assumption of constant Nusselt 
numbers (Nu) seems more appropriate, where: 


Nu = hL 
K 


(15) 


The ratio between the inflow and the conden- 
sation rate must remain constant at corres- 
ponding times in the fill process for therm- 
odynamic similarity to be maintained, hence 


M i(1 « M, 


cond 


(16) 


The characteristic length is proportional to 
the tank radius so: 


^“Widel ^ kfull SC a1e 


(17) 


For the same fluid the thermal conductivity K 
is constant. So: 

N U model ~ Nu fu11 scale (18) 


Implies : 


h 


model 


1 h 
X 


full scale 


(19) 


Seal i ng 

Work on Orbit Transfer Vehicles has 
resulted in the design of tankage for hydrogen 
with volumes of 425 m 3 and required inflow 
rates of 900 kg/hr. To make the testing of 
the 4.96 m 3 RPM tank applicable to OTV design, 
some scaling rationale must be used in select- 
ing test inflow rates. 

A geometrical scale factor can be defined 
assuming tank shape is the same for the OTV 
and RPM tank. 


A inf should be a function of tank geometry only 
(in the same flow regime) so: 



^inf model 

^ ^ inf full scale (20) 

Substituting these 
for the model and 
through 

relations into equation 10 
the prototype and dividing 


~ in, model 

_ kiWiel - \ (21) 


^in.fullscale 

^^ful Iscale 


X 


— ^model— 
R ful1 scale 


3 ^mode 1 
'I V full scale 


( 12 ) 


using equations 14 and 21 yield a test tank 
mass flow rate between 217 kg/hr and 444 
kg/hr. 
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S pray System Selection 

Two flow systems were selected for test- 
ing at K-Site. It was the objective of the 
spray system to bound the No-Vent Fill problem 
by selecting two spray systems which represent 
the best and the worst conditions in zero 
gravity. Current concepts (Chato, 1988; 
DeFelice and Aydelott, 1987; Merino et al, 
1978) of space No-Vent Fill systems use one or 
more pressure atomizing spray nozzles to 
inject the liquid inflow as a stream of drop- 
lets through the ullage, hence, promoting 
condensation of the ullage gas on the droplet 
stream. As the tank fills, these nozzles will 
submerge and it is expected that the outflow 
will transform to a liquid jet within the bulk 
tank liquid. This jet will continue to pro- 
mote condensation by using fluid mixing to 
transport colder liquid to the tank free 
surface. It is expected that the droplet 
spray will produce much higher condensation 
rates than the submerged jet due to the much 
higher surface area available for heat trans- 
fer. 


Unfortunately, for the spacecraft 
designer, the location of the ullage bubble in 
zero-gravity is highly uncertain so prediction 
of the condition under which the spray nozzles 
will submerge is difficult. The selected 
systems for test bound the problem as follows. 
One is a spray system with a single spray 
nozzle at the bottom of the tank. This 
represents the worst case since it will flood 
soon after liquid begins to accumulate in the 
tank. The other spray system uses an array of 
13 spray nozzles spraying from the top of the 
tank (13 was selected due to the availability 
of a commercial spray manifold with this 
configuration) . These nozzles are located in 
a position such that the spray nozzles are not 
submerged even when the tank is 95% full of 
hydrogen (95% is the target fill level for OTV 
operation) . The flow capacities of each 
system are sized so that (as closely as pos- 
sible within the constraints of commercially 
available sizes) the two different flow sys- 
tems have the same inflow rate for the same 
system pressure. Figure 2 shows the position 
of each system in the Tank. 

The test rig is designed to operate as a 
blow down system with an average delta pres- 
sure of 68 KPa . The flow capacity of each 
system was sized to provide near the higher 
required inflow of 444 kg/hr liquid hydrogen 
at a 68 kPa pressure drop. The system is 
designed to operate as a pressurized transfer 
at a constant transfer line pressure of 207 
KPa with tank pressures from vacuum to near 
line pressure so variable flow rates from 45 
kg/hr to 771 kg/hr must be accommodated. Flow 
capacities for commercial nozzles are normally 
given in gallons per minute of water so con- 
version of the required hydrogen flow rate is 
necessary . 


The continuity equation for nozzles can 
be used to convert a mass flow rate into a re- 
quired flow velocity by the following equa- 
tion; 


v = _M (23) 

P A 

For an incompressible liquid the flow velocity 
through a nozzle can be calculated with the 
following : 

V = C D J2II (24) 

Dombrowski and Wolf son (1972) indicate that 

for pressure atomizing nozzles C D is fairly 
constant over a wide range of flow conditions 
for a specific nozzle so C D will be assumed to 
be constant with flow rate and the same whe- 
ther the flow is hydrogen or water. The 

volumetric flow rate is defined as: 

G = Av (25) 

Using equations 23 and 25 the required G for 
hydrogen is 


G H2 = 44 4 kq/hr = 6.2 6 m 3 /hr = 27.6 gpm (26) 
70.8 kg/tir 

Using equations 24 and 26 for the same nozzle 
pressure drop and size 


—M2 = r^water = 3. 76 (27) 

G water J ^H2 

So the design water flow rate at 68 KPa is 
1.67 m/hr (7.34 gpm). 

For the bottom spray a commercial full 
cone nozzle with a flow capacity of 1.89 m 3 /hr 
(8.3 gpm) water at 68 KPa was selected (Spray 
Systems Catalog) . This nozzle has a nominal 
orifice diameter of .9525 cm and a C p of ap- 
proximately 0.6. It is shown in Figure 3. 
For the top spray a manifold of 13 full cone 
nozzles each with a flow capacity of 0.114 
m /hr (0.50 gpm) water at 68 KPa giving a total 
flow of 1.48 m 3 /hr were selected (Spray Systems 
Catalog) . These have a nominal orifice dia- 
meter of 0.20828 cm and a C D of approximately 
0.8. The manifold is shown in Figure 4. 

ANALYTICAL MODELS 

Heat Transfer Model for the Top S pra y 

To model the heat transfer from the 
droplets of the top spray, the correlation of 
Drown (1951) is used. Drown gives: 


hD = 0.667 D 2 
k 4at | 


6 £ 1 e 

“ n=l 772 

n n 


- (4n^7i^at/0 2 ) 


6 D 2 [V - | 1 

47rat 90 n ' 1 n 4 


-(4n 2 7i ? ut/D 2 ) 


(28) 
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The spray is approximated as a monodisperse 
spray with a droplet diameter the same as the 
volume mean diameter of the spray. The manuf- 
acturers mean diameter information for water 
is corrected to hydrogen with the correlation 
of Steinmeyer (1973) ? 

—system ^ fesystem-1 |"^system-l plater 1 ( 29 ) 

^water “ L°water J L^water J L^system J 


For hydrogen 

D H2 = 0.259 (30) 

D wa ^gp 

The droplet residence time is calculated 
as follows: 

t - L/v (31) 


L is estimated as the centerline distance from 
the spray nozzle to the tank free surface. A 
maximum L is the distance from the spray 
nozzle to the tank wall. The maximum L (mea- 
sured off a scale drawing of the tank) was 
combined with a tank fill geometry function to 
give the L values as the tank fills. These 
functions are listed in Table 1. 

Table 1 Top Spray Length Functions 

Top 6 nozzles L = 22.86 cm 

Middle 6 nozzles L = 2x(i58 cm - x),115 cm max 

Centerline Nozzle L = 156 cm - x 

where x is the liquid fill height 

The heat transfer area is determined by cal- 
culating the number of drops in residence from 
the volumetric inflow rate, mean drop diamet- 
er, and flight time. 

n = 6Gt (32) 

ttD 3 

By multiplying by the surface area per drop; 

A jnf = rnrV (33) 

4 

Hea t Transfer from the Bottom Jet 

Dominik (1984) proposes for one-g transfers 
the following correlation for submerged jets 
based on rates established for jets impinging 
on a flat plate. 


hQ = 0.205 (Re) 0 731 Pr 1/3 (34) 

I< 

Where the reference diameter D is the diameter 
of the jet at the free surface. The angle of 
spread used to calculate this diameter is 
based on the work of Idelchik (1986), who 
gives spread angle as 15°. The velocity is also 
measured at the free surface. This is 


estimated as the volume mean jet velocity and 
is given as; 

v — 0.46 (35) 

v 0 0.16 + 0.29 

D o 

The free surface area is given by the fill 
height x and a function of tank geometry as 

A jnf = Jra 2 ( 1 - ( x . - _b)j ) (36) 


similarly fill volume as a function of fill 
height for the RPM tanks is 

V i • = n 2 a 2 x - ira 2 h - rra 2 ( x - b ) 3 (37) 

3 ~3t7 

Where a = 111.25 cm and b = 92.71 cm for the 
RPM tank. In the actual computer algorithm 
this equation is iterated on to find the 
correct fill height given a known fill volume. 

RESULTS 

The NVFILL computer code was modified to 
include the heat transfer and area functions 
given in the previous sections. The solution 
algorithm is the same as original, except the 
equations for heat transfer coefficient and 
interface area (equations 28 and 33 for the 
top spray; equations 34 and 36 for the bottom 
jet) were used to generate values for these 
parameters which were then used to update 
equation 10 continuously. The code was run to 
simulate tests proposed for the RPM tank. 

Figure 5 shows a comparison of pressuri- 
zation rates of top spray fills at different 
initial tank wall temperatures. The first 
temperature 22.2 K is sufficiently low that 
almost no energy remains in the tank wall. 
The second temperature 77.8 K is the same as 
the vacuum chamber cold wall (hence easily 
obtained) . The third temperature 100 K was 
picked to be about as warm a wall temperature 
one could select and still expect to fill the 
tank. As expected the pressure rises quicker 
at the higher wall temperatures. All three 
wall temperatures show acceptable pressure 
rise rates. 

Figure 6 shows a comparison between the 
previous fill with a 22.2 K initial wall 
temperature and a bottom jet fill with the 
same initial wall temperature. The bottom 
fill is much less effective and must be ter- 
minated at only 80% full. Initially the 
bottom jet does a better job of filling by 
preventing stratification. As the fill pro- 
ceeds, the heat transfer rate and available 
surface area decreases causing rapid pressure 
rise. The increasing pressure causes the flow 
rate to drop, decreasing the heat transfer 
rate, and hence even further accelerating the 
rate of pressure rise. Bottom jet fills at 
higher wall temperature were not investigated 
due to the poor performance of the bottom jet 
fill at even the lowest initial wall temper- 
ature. 
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Several analyses were run at a lower 
supply pressure in the expectation that the 
reduced inflow rates obtained would improve 
the no-vent fill process. Figure 7 shows the 
effect of supply pressure on the fill process 
for the top spray system (with an initial wall 
temperature of 22.2 K) . Figure 8 shows the 
same comparison for the bottom fill. Surpris- 
ingly the reduced pressure did not have much 
effect on the final fill state. 

CONCLUDING REMARKS 

A series of tests for a 2.2 m diameter, 
4.96 m 3 liquid hydrogen tank has been proposed 
to study the No-Vent Fill process. Analytical 
models to estimate No-Vent Fill performance 
have been derived for two spray systems. One 
system uses an array of 13 nozzles spraying 
through the ullage space. One system uses a 
single nozzle spraying up from the tank bottom 
which quickly floods producing a submerged 
jet. The veracity of these models will be 
demonstrated by a series of experimental tests 
planned for the summer of 1989. The results 
of experiments will indicate strengths and 
weaknesses of the proposed models and su- 
ggested areas for refinement. The continuing 
refinement of these models will lead to models 
which can be verified by the minimum amount of 
in-space testing which will enable the ana- 
lysis of in-space operational systems. 

Some areas currently under consideration 
for refinement include improved modelling of 
stratification, better simulation of bulk 
boiling, improved gas side heat transfer, and 
the addition of external heat leak effects. 
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FIGURE 1. - 4.96 CU M TANK IN SUPPORT STRUCTURE 
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FIGURE 3. - BOTTOM JET NOZZLE 
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FIGURE *1.-13 NOZZLE TOP SPRAY 
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PRESSURE, kPa PRESSURE, kPa 



FIGURE 5.-13 NOZZLE TOP SPRAY PREDICTED PERFORMANCE. 
EFFECT OF INITIAL TANK WALL TEMPERATURE SUPPLY PRES- 
SURE, 207 kPa; INFLOW TEMPERATURE, 20.3 K, TANK VOLUME 
4.96 Cli M; MASS TO VOLUME RATIO, 30 KG/CU M. 



FIGURE 6. - COMPARISON OF INJECTION TECHNIQUES. SUPPLY 
PRESSURE = 207 KPA; INITIAL TEMPERATURE = 22.2 K, IN- 
FLOW TEMPERATURE, 20.3 K, TANK VOLUME, 4.96 cu M; MASS 
TO VOLUME RATIO, 30 KG/CU M. 
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FIGURE 7.-13 NOZZLE TOP SPRAY PREDICTED PERFORMANCE. 
EFFECT OF SUPPLY PRESSURE, INITIAL TEMPERATURE, 22.2 K; 
INFLOW TEMPERATURE, 20.3 K, TANK VOLUME, A. 96 CU M; 
MASS TO VOLUME RATIO, 30 KG/CU M. 



FIGURE 8. - BOTTOM JET PREDICTED PERFORMANCE. EFFECT OF 
SUPPLY PRESSURE, INITIAL TEMPERATURE 22.2 K; INFLOW 
TEMPERATURE, 20.3 K, TANK VOLUME, 4.96 CU M; MASS TO 
VOLUME RATIO, 30 KG/CU M. 
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